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Scale  ratio  X 


SUMMARY 


The  maximum  force  of  an  ice  sheet  acting  on  a structure  with  vertical  faces  can 
be  written  (for  notations  see  List  of  Symbols)  as 


rmax=rced  <kN> 

assuming  d/e  to  be  so  large  that  the  indentation  effect  is  negligible  (see  Structures 
with  Vertical  Forces). 

The  actual  intermittent  force,  a peak  load  of  short  duration  acting  on  the 
inclined  wedge,  will  be  reduced  to 

f 

CF  is  a function  of  a system  of  nondimensional  data: 


CF 


= 0 


/ £ L X. 
\pu2'd’e 


An  approximate  formula  of  CF,  valid  in  the  interval  0 <;  CF  < 0.4,  can  be  written 
as 


CF  ~ 


s.2>yi~jrc 

Vc 


(A) 


where 


C = 0.1 6 • \/EI(puc2  ■ sin2 a)  ■ C,  /C2  • (C5/C2)2  . 

Equation  A is  derived  by  use  of  a rupture  stress  analysis  and  by  use  of  dynamic 
equations  of  the  system.  It  is  assumed  that  no  other  exterior  forces  are  acting  on 
the  ice  sheet  and  that  uc  > 0. 

The  effect  of  friction  in  the  vertical  direction  on  the  face  of  the  wedge  has  been 
included,  while  friction  in  the  horizontal  direction  has  been  neglected,  as  its  in- 
fluence on  rupture  is  negligible. 

In  the  case  of  uc  = 0,  representing  a floe  resting  against  the  wedge,  the  formula 
becomes 


CF  = 


Ftx  will  be  either  the  wind  shear  or  the  shear  from  a current  or  a combination  of 
the  two,  represented  by  the  formulae,  respectively,  where  A is  the  area  of  the  floe 
(m2): 


^x.w,nd  = 4.8-10-3-ipa  t;2,r,d  /1  (kN) 

/L.x.«lnd  = S.4-10-3-lpvv-/:/2ater-4  (kN) 

where  pa  = 1 .25  x 1 0'3  tonne/m3  is  the  density  of  air  and  pw  = 1 tonne/m3  is  the 
density  of  water,  which  gives  a value  of  CF  = 1 .0  for  /~ex  = rced.  For  Fex  > rced 
the  floe  will  move  continuously  and  eq  A will  apply. 

The  velocity  range  will  be  from  approximately  0.1  m/s  to  4 m/s.  The  angle  a 
may  vary  from  30°  to  60°  and  (3  from  45°  to  70°.  The  friction  coefficient  will 
seldom  exceed  0.1 , but  in  some  cases  it  may  be  considerably  higher.  The  value  of 
C,  /C2  will  vary  from  0.1  to  1 .0  depending  on  the  combination  of  a,  P and  p. 
Graphs  showing  the  relationships  are  presented  in  Figures  7 and  10.  The  factor 
e/d  can  vary  from  zero  to  approximately  0.3,  limited  by  the  assumptions  made 
for  the  rupture  patterns.  The  ratio  between  the  bending  strength  and  the  com- 
pression strength  may  vary  from  0.2  to  0.5. 

The  influence  of  p,  uc  and  E is  not  appreciable,  since  it  is  in  the  fourth  root. 
However,  an  increase  in  velocity  should  give  an  increase  in  force.  Low  values  of 
E should  also  result  in  higher  values  of  CF.  An  increase  in  P will  result  in  an  in- 
creased force.  If  p -*■  0.3,  C^/C?  can  reach  rather  low  values,  resulting  in  values  of 
CF  up  to  1 .0,  i.e.  no  reduction  at  all.  The  theoretical  curves  and  the  results  of  the 
tests  are  shown  in  Figure  1 2. 
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INTERMITTENT  ICE  FORCES  ACTING 
ON  INCLINED  WEDGES 

Per  Tryde 

INTRODUCTION 

The  research  described  in  this  report  was  initiated  at  the  Institute  of  Hydrodynamics  and 
Hydraulic  Engineering  (ISVA),*  Technical  University  of  Denmark,  by  Professor  H.  Lundgren,  as 
methods  for  the  rational  calculation  of  ice  forces  acting  on  structures  in  the  open  sea  were  urgently 
needed. 

In  1971  a simple  formula  for  calculating  ice  forces  was  developed  for  structures  with  vertical 
faces  which  were  either  perpendicular  to  the  direction  of  the  forces  or  positioned  as  wedges  with 
vertical  faces.  It  was  soon  realized  that  the  design  forces  could  be  considerably  reduced  if  inclined 
faces  were  introduced,  thus  resulting  in  bending  failure  instead  of  compression  or  shear  failure. 
Korzhavin  (1971)  presented  an  empirical  formula  but  did  not  establish  the  functional  relationship 
of  the  parameters  involved  in  the  system.  The  aim,  therefore,  was  1)  to  present  a theory  that  could 
explain  the  phenomenon  and  2)  to  verify  the  theory  by  model  tests. 

The  theory  was  presented  in  1973,  and  has  since  been  improved  in  accordance  with  the  findings 
from  model  tests.  Model  tests  were  made  in  Lyngby  with  an  artificial  "ice”  material  developed  by 
ISVA.  Tests  with  natural  ice  have  been  performed  by  the  Hamburgische  Schiffbau-Versuchsanstalt 
(HSVA)  in  Hamburg,  Germany,  and  by  CRREL  in  Hanover,  New  Hampshire. 

RECENT  ICE  RESEARCH  IN  DENMARK 
General 

This  section  contains  a short  review  of  pertinent  ice  research  performed  at  ISVA  since  1968. 

In  Denmark,  ice  forces  had  hitherto  been  determined  by  empirical  rules  based  or:  the  Danish 
Standards  of  Loading,  indicating  a horizontal  ice  force  on  a bridge  pier  corresponding  to  1 .5  to 
3.0  tonnes  per  meter  free  span  of  the  bridge  in  question,  the  higher  value  to  be  used  for  structures 
in  fiords,  sounds,  etc. 

Recently,  construction  activity  in  Danish  waters  has  increased  with  the  building  of  large  bridges 
interconnecting  the  main  islands,  oil  piers  on  the  open  shores,  and  offshore  lighthouses. 

Structures  with  vertical  faces 

A general  expression  for  ice  force  can  be  written  as 

F = krced 

* Instituttet  for  Stromnlngsmekonlk  og  Vandbygnlng. 


where  k is  a dimensionless  variable,  rc  the  compressive  strength  of  the  ice,  e the  ice  thickness,  and 
d the  width  of  the  structure.  For  slender  vertical  structures  such  as  piles  and  cylinders  the  value  of 
k can  be  given  as 


l 

r 


A = — = 1 + 1.5%  for  0 «i  % 2 

o„  d d 


where  o„  is  the  strength  for  the  case  d -*  This  strength  may  not  be  equal  to  what  is  normally 
denoted  as  the  compressive  strength,  since  firm  rules  for  performing  “compression"  tests  are  lacking. 
This  problem  is  at  present  being  considered  by  the  IAHR  committee  on  ice  testing.  The  formula 
for  the  force  F may  therefore  have  to  be  corrected  accordingly.  The  effect  of  ice  thickness  on  com- 
pressive strength  has  been  studied  by  Schwarz  (1974),  and  by  Frederking  and  Gold  (1974). 

The  formula  k - 1 +1.5  e/d  will  give  values  which  are  too  high  when  e/d  is  larger  than  1 , 
according  to  investigations  by  Afanas'ev  (1973),  Assur  (1972)  and  Schwarz  (1974).  The  following 
formula  will  give  results  between  those  of  Assur  (1972)  and  Afanas’ev  (1973): 


* = 1+2.1 


-1 


Wedges  with  inclined  faces 

Korzhavin  (1971)  has  given  formulae  for  ice  forces  acting  on  inclined  wedges,  based  on  analytical 
considerations  and  data  obtained  from  prototype  structures.  However,  a theory  has  not  been  pre- 
sented by  Korzhavin.  The  author  attempted  to  develop  such  a theory  and  presented  this  in  a lecture 
given  at  the  University  of  Iowa  Institute  of  Hydraulic  Research  in  the  summer  of  1973,  the  sub- 
stance of  which  is  described  in  Tryde  (1973b). 

To  evaluate  this  theory,  model  tests  were  performed  at  ISVA.  These  tests  were  made  in  a 2-m 
wide  flume  with  a floe  of  approximately  8 m2,  of  which  approximately  2.5  m2  consisted  of  an 
artificial  material  composed  as  stated  below  and  the  remaining  part  consisted  of  plywood.  The 
model  floe  was  run  against  a wedge-faced  object  on  a flowing  current.  In  these  tests  the  Froude 
number  was  approximately  0.2,  with  a floe  velocity  of  0.3  m/s,  an  ice  thickness  of  approximately 
1 cm,  and  a wedge  width  of  16  cm.  The  object  had  an  included  angle  at  the  point  of  the  wedge  of 
100°  in  the  horizontal  plane  and  an  inclination  of  55°  to  horizontal.  The  wedge  could  be  rotated 
± 10°  in  order  to  vary  the  angle  of  inclination.  The  wedge  was  mounted  on  a dynamometer  equipped 
with  strain  gauges,  making  it  possible  to  record  the  forces  acting  on  it. 


Recording  of  the  strength  of  natural  ice 

In  subsequent  tests  it  is  planned  to  perform  in-situ  measurements.  Bending  and  compression 
tests  will  most  likely  be  selected  as  a relevant  determination  of  the  strength  parameters.  A loading 
frame  has  been  developed  in  order  to  test  beams  while  they  are  still  afloat  in-situ. 


ARTIFICIAL  ICE  - PROPERTIES  AND  MANUFACTURE 
General 

In  connection  with  research  concerning  ice  forces  acting  on  structures,  it  was  necessary  to  perform 
model  tests  in  order  to  verify  various  theories.  Since  cooling  facilities  were  not  available  at  ISVA  it 
became  essential  to  develop  an  artificial  material  with  properties,  at  temperatures  around  20°C,  con- 
forming to  the  requirements  of  the  relevant  model  laws. 
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Figure  I.  Grading  diagram  of  granu-  Figure  2.  Setting  time  of  plaster  of  Paris  for  van- 

late  of  plastic  described  in  this  report.  ous  salt  solutions. 


Investigations  of  floes  drifting  against  structures  involve  dynamic  actions,  leading  to  the  Froude 
model  law  comprising  similitude  of  inertia  forces  of  water  and  ice,  internal  forces,  and  gravity  forces. 
[Froude  number  = ujs/g  ■ e,  where  g is  acceleration  of  gravity  (m/s2),  and  uc  is  the  velocity  (m/s) 


of  the  floe.) 

The  scale  of  the  models  was  as  follows: 

Geometric  (model  to  prototype)  X = 1m/1 

Velocity  \/\ 

Acceleration  1 

Time  v/X 

Force  X3 

Mass  X3 

Density  1 

Pressure  X 

Strength  X 

Young’s  modulus  X 

Friction  1 . 


With  model  scales  of  the  order  1/10  to  1/50  this  led  to  a material  with  a bending  strength  of  10 
to  100  kN/m2,  corresponding  to  a bending  strength  of  100  to  5000  kN/m2  for  ice  in  nature. 


Composition  of  material 

The  material  was  composed  as  follows: 

Plaster  of  Paris  (CaS04)  600  g 

Granulate  of  plastic,  1-2  mm  diam  (light  weight)  80  g 

Salt  (kitchen  salt)  100  g 

Borax  (Na2B407)  2g 

Air-entraining  agent  10  g 

Water  300  g 

which  yields  approximately  one  liter  of  material. 


It  should  be  noted  that  the  above  quantities  may  require  adjustment  to  compensate  for  variations 
in  the  properties  of  the  plaster  of  Paris.  The  plastic  granulate  was  expanded  by  the  Danish  firm 
ANWI  from  polystyrene  supplied  by  Sinclair-Kopper  of  Pittsburgh,  Pennsylvania.  The  granulate  has 
a grain  diameter  between  0.5  and  2.0  mm  for  ice  thicknesses  between  1 and  2 cm  (see  Fig.  1 ).  This 
can  be  varied  as  required. 

The  unit  weight  of  the  granulate  is  from  40  g/liter  up  to  approximately  80  g/liter.  It  is  possible 
to  re-use  the  granulate  by  crushing  the  sheets  and  floating  the  plastic  from  the  plaster  of  Paris  paste. 
There  will  be  a slight  increase  in  unit  weight  by  this  process,  which  may  be  taken  into  account  when 
composing  the  material. 
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Salt  is  used  to  reduce  the  strength  of  the  material.  A saturated  salt  solution  produces  very  low 
strength.  Borax  acts  as  a retarder,  delaying  the  setting;  a 0.3%  solution  gives  a setting  time  of  about 
I S min  (at  20°C).  Eight  minutes  is  the  normal  setting  time  for  plaster  of  Paris.  Borax  will  counter- 
act the  accelerating  effect  produced  by  the  salt  (see  Fig.  2).  The  air-entraining  agent  used  is  a resin 
giving  approximately  4%  air  content. 

The  hardened  product  saturated  with  water  has  a density  of  0.90  to  0.95  g/cm3  and  a bending 
strength  of  70  to  200  kN/m2.  The  compressive  strength  is  3 to  5 times  as  great.  Young’s  modulus 
is  from  105  kN/m2  to  106  kN/m2,  and  friction  p = 0.1  to  0.3  against  a polished  plastic  plate. 

Order  of  mixing  ingredients 

The  ingredients  of  the  artificial  ice  material  are  mixed  in  the  following  order: 

1 . Make  the  salt  solution 

2.  Add  the  air-entraining  agent 

3.  Add  plaster  of  Paris 

4.  Add  granulate  of  plastic. 

A small  part  of  the  salt  solution  should  be  retained  for  final  adjustment  of  the  plasticity  of  the 
mix. 


RUPTURE  THEORY 

In  some  problems  involving  the  destruction  of  materials  it  is  important  to  study  the  type  of 
rupture,  which  is  visible  in  the  form  of  rupture  lines  and/or  zones.  The  theory  of  elasticity  possesses 
the  inherent  defect  of  not  giving  any  information  about  safety  against  ultimate  failure,  although 
in  practical  engineering  there  seems  to  be  an  increasing  tendency  to  employ  design  methods  which 
are  based  on  the  state  of  failure.  Coulomb’s  method  was  probably  one  of  the  first  to  use  this 
principle.  However,  some  problems  are  too  complicated  to  be  solved  by  use  of  the  theories  of 
elasticity  unless  very  simplifying  assumptions  are  made,  which  sometimes  cause  inaccurate  results. 

By  using  the  theory  of  plasticity  it  is  possible  to  determine  the  lines  and  zones  of  rupture  - in 
principle.  In  practice  this  may  present  some  other  problems  as  exact  integration  can  only  be  carried 
out  in  a few  simple  cases.  However,  this  can  sometimes  be  overcome  by  the  use  of  computers. 

Dynamic  action  introduces  further  complications  in  establishing  the  location  of  rupture  lines  or 
zones.  Therefore,  in  the  present  investigation,  elastic  as  well  as  plastic  theory  will  be  used  in  order 
to  determine  the  unknown  parameters  of  the  rupture.  This  is  most  easily  demonstrated  as  the  dif- 
ferent formulae  are  derived. 

The  basic  ideas  of  this  investigation  will  first  be  tried  on  a simplified  case  of  an  tnclined  plane, 
where  solutions  can  be  obtained  by  use  of  the  theory  of  elasticity.  Then  they  will  be  compared 
with  the  method  based  on  plasticity  and  with  empirical  methods  (Reeh  1972,  Tryde  1972  and 
Sorensen  1976). 


THEORETICAL  APPROACH 
Initial  contact 

When  a floe  drifts  against  an  inclined  wedge  (see  Fig.  3)  initial  contact  is  made  at  the  ridge,  where 
the  lower  edge  of  the  ice  sheet  is  crushed,  while  the  floe  is  pushed  along  the  ridge. 

A force  acts  at  the  point  of  contact  (position  I in  Fig.  3)  and  the  vertical  component  of  this  force 
will  produce  a crack  propagating  along  the  main  axis  of  the  wedge.  The  force  necessary  to  produce 
this  initial  crack  can  be  calculated  from  the  theory  of  an  elastic  plate  on  an  elastic  support.  As  the 
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Figure  3.  Plan  and  section  of  wedge  ( where  a is  one-half 
the  angle  at  the  point  of  the  wedge  in  the  horizontal 
plane,  V the  vertical  force  and  F the  horizontal  force, 
both  without  friction). 


Figure  4.  Location  of  rupture  in  front  of  wedge. 
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floe  moves,  contact  is  established  along  the  faces  of  the  wedge,  and  it  can  be  shown  that  the  area 
adjacent  to  the  wedge  will  be  a rupture  zone.  This  indicates  that  rupture  will  appear  randomly, 
either  as  crushing  along  the  face  or  as  a combined  compression/bending  failure  at  some  distance 
from  - and  parallel  to  - the  face.  During  this  phase  an  increasing  area  of  the  floe  will  come  in 
contact  with  the  wedge  and  the  force  will  gradually  increase  until  the  floe  engages  the  full  width  of 
the  wedge. 

By  the  time  the  floe  has  reached  position  II  in  Figure  3 (ensuring  complete  penetration),  the 
cracking  pattern  will  adjust  itself  to  a stage  of  repeated  occurrence.  The  crack  in  the  centerline 
will  propagate  continuously,  while  the  parallelograms  in  front  of  the  wedge  will  break  off  at  a given 
distance  from  the  face  (see  Fig.  4). 

The  initial  vertical  force  P , required  to  produce  the  crack  can  be  found  using  the  theory  given 
by  Nevel  (1965).  Provided  the  force  is  acting  on  a small  area  of  the  ice  along  the  ridge: 


*r 


a 1.0 


where  Mx  is  the  unit  moment  along  the  crack  at  the  edge  and  P , is  the  It,-.  J applied  at  the  edge. 

The  rupture  moment  can  be  written  = rbW,  where  rb  is  the  b eng  strength  and  H'  is  the 

section  modulus.  The  magnitude  of  the  force  will  thus  be 


where  e is  the  thickness  of  the  ice.  As  an  example,  rb  = 500  kN/m2  and  e = 0.5  m gives  Py  = 21  kN, 
which  is  the  initial  vertical  force  necessary  to  rupture  the  ice. 
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Figure  5.  Initial  stage  of  rupture 
along  edge  of  floe. 
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Vertical  Section  2-2  Section  3-3 


Figure  6.  Plan  and  sections  of  inclined  wedge  and  act- 
ing forces  (dotted  lines  indicate  case  for  p = 0). 


After  the  formation  of  the  longitudinal  crack,  the  ice  along  the  faces  will,  in  the  initial  stage,  be 
crushed  along  the  faces  of  the  wedge.  This  is  illustrated  in  Figure  5 where  'A  P2  is  the  upward 
vertical  force  required  to  break  off  each  corner.  The  section  modulus  along  A-B  will  be  2 Z e2/ 6 
where  7.  is  the  distance  from  the  point  of  contact  out  to  the  crack.  This  will  give  P2  = rb  e2  2/3, 
i.e.  a force  which  is  four  times  as  great  as  /*, , both  sides  breaking  off  simultaneously. 

This  situation  only  prevails  during  a short  period  before  the  parallelogrammic  pattern  is  initiated. 
The  force  is  independent  of  the  width  of  the  wedge,  and  for  all  practical  purposes  it  is  rather  small, 
being  only  a small  percentage  of  the  total  force  for  which  the  structure  must  be  designed. 


Forces  acting  on  wedge 

The  geometry  and  the  principal  force  diagrams  for  forces  acting  on  a wedge  are  illustrated  in 
Figure  6.  In  this  figure  N is  the  force  normal  to  the  face,  T’  is  the  force  in  the  horizontal  plane, 
is  the  vertical  force  on  the  face  and  0 is  the  inclination  of  the  wedge  to  horizontal.  The  coef- 
ficient of  friction  is  denoted  as  p,  and  all  forces  with  a prime  mark  indicate  that  friction  has  been 
included.  The  forces  can  be  expressed  as  follows: 

r*1'- 


where 


= uni 

sina 


F=2T‘  sina  = 21 f,C2  sina 

v|  = V,  jjjo£. ) - V,  c, 


where 
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c, 


= 1 


-uMli 

sina 


and  Vy  is  the  vertical  force  on  each  face  in  case  of  no  friction.  The  functional  relationship  of  the 
various  parameters  is  shown  in  Figure  7. 

Friction  in  the  horizontal  plane  has  been  neglected,  as  the  influence  of  this  on  the  rupture  is 
negligible.  The  friction  in  the  vertical  direction  may  possibly  be  increased  by  ice  being  frozen  to 
the  face  of  the  wedge  (thus  preventing  the  ice  from  moving  upwards).  Flowever,  as  the  rupture 
patterns  will  hardly  be  affected,  the  increase  of  the  total  force  will,  for  p = 0,1 , amount  to  approxi- 
mately 10%. 

Maximum  force  and  actual  force 

The  maximum  force  Fmax  on  a given  section  in  front  of  the  structure  is 


where  rc  = compression  strength  (kN/m2) 
e = thickness  of  the  ice  (m) 
d = width  of  the  structure  (m). 
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Figure  8.  Rupture  lines  in  front  of  inclined  wedge  (b  = width  of 
rupture  channel). 


The  actual  force  is  expressed  as 


F-CrF^ 


where  CF  is  a reduction  factor,  which  is  a function  of  a system  of  nondimensional  data: 


where  £ = Young's  modulus  (kN/m2) 
p - density  of  ice  (tonnes/m3) 
uc  - velocity  of  the  floe  (m/s) 
y - characteristic  breakoff  distance  (m) 
rb  = flexural  strength  (kN/m2) 

a = included  angle  at  the  point  of  the  wedge  in  the  horizontal  plane 
0 = inclination  of  wedge  to  horizontal 
p = coefficient  of  friction. 

Principal  stress  equations 

The  assumed  rupture  lines  are  indicated  in  Figure  8. 

The  load  along  the  face  on  each  side  of  the  wedge  is  assumed  to  be  uniformly  distributed.  The 
vertical  and  the  horizontal  forces  produce  an  eccentric  load  combination  on  the  rupture  plane  in  the 
ice,  which  results  in  the  following  maximum  boundary  compressive  stress: 


v\y  - T’x 


T'(kyyl  Una) 


e [bl 2 sin  a)  i/4  e7  {b/2  sina)  >/*  [b/2  sin  a)2  e 

where  y is  the  breakoff  distance  as  shown  in  Figure  8,  and  x is  the  eccentricity  of  the  horizontal 
force,  which  can  be  expressed  as  x = fe,  f being  the  eccentricity  coefficient.  This  coefficient  can  be 
written  as  f = Vi  (1  - Cr ).  The  eccentricity  in  the  horizontal  plane  would  be y/ tana  if  the  width 
of  the  rupture  channel  were  equal  to  the  width  of  the  wedge.  As  additional  cracks  develop  as  shown 
in  Figure  8,  it  is  proposed  to  adjust  the  eccentricity  to  a value  between  the  two  extremes  illustrated 
in  Figure  8.  Therefore,  the  factor  h1  has  been  introduced. 
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For  the  section  modulus  'It  (bl2  sina)2  e,  a modified  expression  1.5  • '/«  \dbe/(2  sinor)2  j will 
be  used  in  order  to  simplify  the  following  formula: 


_ F,  Cj  2 sina  ^ (V,  C,y  - F',  C2x)  12  sina  (/  r k,y  24  sin2a 
0bc  eb  eJfj  1 2 T.Sbde  Una 


(2) 


using  the  forces  derived  in  the  section  Forces  Acting  on  Wedge. 
By  inserting  x = (e/2)(1  - CF ),  eq  2 can  be  written  as 


be 


2 K,  sina 
eb 


C2  |1  - 3(1  -CF) | 6 C,  +2L  C2  8b,  -L  cosa 


(3) 


Introducing  the  system  constant  of  C5  = 6 C,  + bC2  (e/d)  cosa  and  a value  oi  k,  - 0,75,  eq  3 can  be 
written  as 


2 V.  sina 

0bc  7b 


C2(3CF-2)+f  C5|  . 


(4) 


By  using  F‘  - 2 V,  C2  sina : 

V,  = F' ? 

1 2 sina  C2 


cf-r 


F'  _ F' 


rbcebC2 


rced  ~ rced\C2(3CF  - 2)  + (y/e)  C5)  ' 


(5) 


The  rupture  strength  r ^ is  a function  of  e (the  rb/rc  ratio)  and  the  reduction  factor  CF  (see  Fig. 
9).  It  can  be  shown  that  the  rupture  strength  for  a combined  axial  load  and  bending  moment  will 
be  a tension  failure  in  the  interval  0 < CF  < (1  - e/2)  giving 


= 2Cf  ♦ e 

r 


(6) 


and  a compression  failure  in  the  interval  1 - e/2  < CF  < 1 .0  giving 


Figure  9.  Rupture  stress  characteristic. 


— = 1 . 


(7) 


For  the  width  b of  the  rupture  channel  the  following 
empirical  formula  is  suggested,  where  the  ratio  b/d 
is  a function  of  CF : 


£ =-  1.5Cf  +2.5 
d 


(8) 


which  is  based  on  experimental  data.  It  is  correct 
for  CF  = 1,  corresponding  to  the  case  with  vertical 
faces.  Equation  5 will  give,  when  inserting  eq  6, 

7 and  8: 
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(9) 


- 3 C%  + SCF  - 1.5eCF  + 1.5e  + e 
3CF  -2+  (y/e)(CB/C2) 


or  as  an  expression  oi y/e 


(f)r 


- 6CF  + 7 - 1 .5  e + (2.5  e/CF) 


for  interval  1 

(o<cF<i^ 


r - - 1 .5  CF  + 2.5 

F 3Cf  - 2 + (y/e)(C6/C2) 

or  as  an  expression  of  y/e 

/x\  _ ~ 3 CF  + 0.5  + (2.5/Cf  ) 

\ e /2  (C5/C2) 


for  interval  2 

(!^!<Cf<1.o) 


Two  unknowns  appear  in  this  equation,  a solution  requiring  that  one  more  equation  be  found.  This 
is  done  further  below,  where  the  dynamic  equation  is  derived  for  the  ice  pieces  breaking  off  in  front 
of  the  wedge. 

Equation  10  may  be  used  to  determine  the  value  of  y/e  in  the  case  of  a floe  at  rest  in  front  of 
the  wedge,  as  we  then  know  the  values  of  the  force  F equal  to  the  wind  shear  or  the  water  shear 
acting  on  the  floe,  and  Fmax  is  also  known.  Thus 

Wind  shear:  t,  = 4.8  • 10'3  • jpa  10-6  L%,lnd  kN/m2 

Water  shear:  rw  = 5.4  • 10"3  10'3  fT2^,^  kN/m2 

where  l/wind  and  are  the  wind  and  water  velocities,  respectively.  For  F = rced  we  obtain 

CF  = 1.  For  F > rced  the  floe  will  move  continuously. 

To  determine  the  values  of  C,/C2  and  C5/C2  = C3  use  Figures  7 and  10. 

Dynamic  equations 

The  equation  of  angular  momentum  of  the  ice  piece  broken  off  can  be  written  as 

l(yV\-r'x)dt  = lAdu>  (12) 

where  /A  is  the  moment  of  inertia  about  A-A,  x the  eccentricity  of  horizontal  force,  and  to  the 
angular  velocity.  The  vertical  force  can  be  written  (see  Fig.  1 1 ) as 

2h  (1 3) 

1 2 sina  h ' 


which  gives 
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V\  2 sina 


The  time  required  for  making  an  indentation  zh  is 


*h  . 21/;  _ 2 VyCy 


uc  sina  uc  rcd  ucrcd 


The  moment  of  inertia  about  A-A  (in  Fig.  1 1 } will  be 


/A  = ca  ]?y3e  yfe-- 

A A 3 2 sina 


Actually  the  moment  of  inertia  is  slightly  larger,  due  to  the  increased  width  b of  the  rupture  channel. 
CA  will  be  explained  below. 

From  the  theory  of  elasticity,  at  time  dt  the  deflection  of  the  ice  sheet  measured  vertically  at 
the  point  of  contact  is 


U=cB 


Vy  (C,y  -C2ef)y7 


where  f = 1/2  (1  - CF).  The  second  moment  of  area  of  section  A-A  is 


/ = _L  e3 


* 12  2 sina 


Inserting  /,  in  eq  17: 


fA  ~ CB 


8 Vy  ( Cyy  - C2ef)y 2 sina 
Ee3b 


CB  will  be  explained  below.  The  effect  of  the  horizontal  force  on  the  deflection  is  not  exactly  as 
stated,  but  the  error  is  small  and  is  justified  by  the  simplification  obtained.  CA  was  introduced  in 
eq  17  to  allow  for  the  effect  of  the  virtual  mass  of  the  water  below  the  ice  (Fngelund  1966). 

The  virtual  mass  effect  may  be  derived  as  follows.  For  a strip  1 m wide  the  virtual  moment  of 


inertia  /„  is 


/ ~0.2py4  (from  Engelund  1966). 


If  it  is  assumed  that  y ~ 3e  the  moment  of  inertia  of  the  strip  of  ice  rotated  would  be  approximately 


ii-T  py3e  = 1 py*  ~~ 01  py* 


lv*t  = lv*li=0.3py4  = 3l,. 


Actually  the  resulting  moment  of  inertia  is  a function  of  y/e,  but  for  simplification  we  shall  use  a 
factor  CA  = 3 for  the  moment  of  inertia  of  the  ice  piece  itself  in  order  to  obtain  the  total  moment 
of  inertia. 


If  full  fixity  is  assumed  at  the  rupture  line  (A-A),  the  value  of  CB  will  be  1 .0.  If,  on  the  other 
hand,  one  assumes  that  the  action  is  that  of  a continuous  beam  with  the  next  support  at  a distance 
4 y from  the  line  of  rupture,  the  deflection  will  be  five  times  larger,  giving  CB  = 5.0,  which  is  a 
realistic  assumption.  The  effect  of  these  values  is  not  as  great  as  it  would  seem  at  first  sight,  as  they 
appear  in  the  1/4  power  in  the  final  formula. 

The  instantaneous  change  of  cu  may  be  expressed  as 


U _ _ 4(C^y -C2ef)yucrcs\nad 

au)  - ~r — = Cn  • 

dty  B Ee*Cyb 


(23) 


Inserting  eq  1 5 and  eq  1 9 is 


l/,  = 0.82  x/C^CfT  yff  ^ 


£ _L  u r £- 


d± 

b ' 


(24) 


Inserting  in 


2 Vt  sina 


be  eb 


C2(3CF-2)+f  C6 


(25) 


ubc 


1.64  x/CACB  \fp[E  (C2/C,)  rcuc  sina 


by  using 


'be  _ 

r. 


= 2Cc  + e 


and  -te-  = 1 
rc 


in  interval  1 ^0  CF  < * ^ e J , 
in  interval  2 ^ < CF  c 1 ^ . 


(26) 


Furthermore,  if  b/d  = - 1 .5  CF  + 2.5  the  following  are  found: 
interval  1 

o<"  fsr 

interval  2 


(27) 


061  V5F 

■ , 1 0c,-2).S«rf)3 

v/(-1.5CF  + 2.5)3  |Wa  C2  \e  / 2 


(28) 
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in  which  eq  10  and  1 1 should  be  inserted.  When  applying  C'A  = 3 and  CB  = 5 in  eq  16  and  19,  a 
system  parameter  C can  be  introduced  for  intervals  1 and  2,  respectively: 


and 


C = °.,6  J L-  (Ji)2 

W pul  sin?a  l2  \l2  / 
2CF+e  \[e),  \c2  / 


V(-  I.SCp  + 2.5)® 


(3Cf  - 2|  ♦ (f  )] 
(3CF-2).(fj“ 


(29) 


(30) 


(31) 


Intermittent  nature  of  force 

The  acting  forces  are  of  an  intermittent  nature  since  they  increase  abruptly  each  time  a 
parallelogrammic  piece  is  about  to  break  off,  after  which  the  forces  decrease  until  the  next  encounter 
takes  place.  The  period  of  time  between  each  peak  can  be  expressed  by 


tc=-Xr—  = (A-^~  (32) 

c uc  sina  \ e J uc  sina 

where  y/e  can  be  taken  from  eq  10  and  1 1 . Each  peak  of  the  force  takes  place  within  1/5  to  1/50 
of  a second. 


Theoretical  curve 

The  curves  illustrating  the  functions  in  the  two  intervals  are  shown  in  Figure  12.  The  results 
of  mode!  tests,  which  will  be  described  later,  have  also  been  plotted  (see  Plotting  of  Results  below). 
The  dispersion  of  results  is  rather  noticeable,  but  the  agreement  is  fairly  good. 
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An  illustrative  approximation  of  these  theoretical  curves  can  be  expressed  as 


for  interval  1 


0 < CV 


(33) 


VERIFICATION  OF  THEORY  BY  MODEL  TESTS 
Description 

In  order  to  verify  the  theoretical  investigation,  model  tests  using  the  artificial  material  described 
in  Artificial  ice  - Properties  and  Manufacture  were  performed. 

A flume  2 m wide  with  a water  depth  of  0.3  m was  used.  The  experimental  floe,  with  an  area 
of  approximately  8 m2,  consisted  of  plywood,  except  for  an  area  of  1 .2x2.0  m,  where  the 
artificial  material  was  used  (sec  Fig.  13). 

An  inclined  wedge  0.16  m wide,  with  a = 50° 
and  0 = 55°,  was  used  as  the  object.  The  whole 
wedge  could  be  rotated  in  order  to  test  various 
values  of  0.  The  wedge  was  made  of  smooth 
Perspex  in  order  to  minimize  friction.  A 
dynamometer  with  strain  gauges  was  placed  in- 
side the  wedge,  making  it  possible  to  record 
instantaneous  horizontal/vertical  forces  and  the 
moment.  These  force  and  moment  signals  as 
well  as  a velocity  signal  in  the  form  of  pulses  from  a photo  cell  beam  intercepted  by  a “comb" 
placed  on  the  floe  were  recorded  on  a high  speed  (30-in. /s)  tape  recorder. 

For  visual  observation  of  the  signal,  the  tape  recorder  was  played  back  at  slow  speed  (1 5/32 
in./s)  and  a 5-channel  pen  recorder  was  used  to  obtain  curves  from  all  recording  components  (see 
Fig.  14). 

Tests  were  made  with  varying  velocities.  For  detailed  information  concerning  test  results,  see 
App.  A. 

Width  of  rupture  channel 

Test  results  show  that  the  width  of  the  rupture  channel  is  considerably  larger  than  the  width  of 
the  wedge,  in  some  cases  about  2.5  times  as  great.  It  has  not  been  possible  to  develop  a theory  that 
can  predict  the  rupture  pattern.  Therefore,  the  width  is  given  by  the  empirical  formula  b/d  = 

-1.5CF  + 2.5,  which  gives  fairly  good  agreement  with  the  tests. 

Plotting  of  results 

The  moment  and  horizontal  and  vertical  forces  were  recorded  during  the  break  off  of  each  of 
the  parallelogrammic  pieces.  The  CF  coefficient  is  determined  from 
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F igure  13.  Plan  of  flume. 


where 


which  is  known  for  each  test,  as  the  compression  strength  and  thickness  of  the  ice  are  measured. 
Furthermore,  the  velocity  of  the  floe  was  recorded  and  the  size  of  the  breakoff  measured. 
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Figure  14.  Typical  force  variation. 

Each  test  comprised  about  20  encounters.  Based  on  the  magnitude  of  the  vertical  and  horizontal 
forces  and  the  geometry  of  the  wedge,  it  was  possible  to  determine  the  friction  coefficient.  The 
variation  of  the  friction  was  usually  not  appreciable;  therefore,  a mean  value  was  adopted.  Thc 
system  data  C,/C2>  C5/C2,  C were  determined  for  each  test.  At  the  first  encounter,  only  where 
the  floe  had  been  precut  with  straight  lines  to  fit  the  wedge  plan  profile  was  it  possible  to  obtain 
a simultaneously  symmetrical  breakoff.  In  no  other  encounter  during  the  tests  had  this  been 
experienced.  It  was  therefore  concluded  that  symmetrical  occurrence  is  rare. 

Depending  on  the  regularity  of  the  rupture  line,  which  is  a straight  line  only  in  theory,  the  rising 
time  of  the  load  is  approximately  1/100  s.  During  that  time  the  floe  drifting  with  a velocity  of,  for 
examplt  50  cm/s,  would  move  0.5  cm,  which  is  the  size  of  the  deviation  from  the  straight  line  of 
the  rupture.  1 herefore,  it  is  likely  that  the  two  peaks,  on  each  side  of  the  wedge,  will  occur  in  a 
staggered  pattern.  This  implies  that  the  recorded  CF  value  from  the  test  should  be  multiplied  by  a 
mean  value  of  1 .3  before  being  plotted.  In  Figure  12,  the  test  results  have  been  plotted  accordingly, 
so  as  to  compare  them  with  the  theory  based  on  simultaneous  breakoff. 

As  shown  in  the  recordings  of  the  forces,  the  rupture  force  shows  a clear  regularity  although  the 
rupture  pattern  itself  may  look  rather  irregular.  The  forces  acting  are  of  an  intermittent  nature, 
since  they  increase  abruptly  each  time  a parallelogrammic  piece  is  about  to  break  off,  after  which 
the  forces  decrease  until  the  next  encounter  takes  place.  Usually,  it  is  observed  that  the  rupture 
is  asymmetrical,  as  the  ice  pieces  on  each  side  of  the  wedge  do  not  break  off  simultaneously  but 
rather  in  rapid  succession  after  each  other  (see  Fig.  14). 

A few  tests  have  been  performed  with  natural  ice  at  Hamburgische  Schiffbau-Versuchsanstalt. 

The  results  are  in  fair  agreement  with  the  theory,  although  the  ruptures  apparently  are  slightly 
different  from  those  observed  with  artificial  ice.  The  natural  ice  was,  however,  weaker  than  the 
artificial  material  used  in  Copenhagen.  It  appears  that  the  parallelogrammic  pieces  were  broken  up 
into  smaller  parts  at  Hamburg,  but  whether  this  took  place  simultaneously  was  difficult  to  observe. 
This  should  be  further  studied. 
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ENGINEERING  APPLICATION 


The  theory  presented  in  the  preceding  pages  is  impractical  for  engineering  calculation,  because 
the  formulae  arc  complicated  and  lack  a clear  presentation  of  the  functional  relationship  of  the 
various  system  constants.  It  is  therefore  recommended  to  use  the  formula: 


The  forces  acting  are,  as  stated,  of  an  intermittent  nature.  The  period  of  time  between  each  peak 
force  can  be  expressed  as 


'c  = - tMt-  yc  -4— 

c uc  sina  \ e ) uc  sina  C5/C2  sina 

where  an  approximate  equation  y/e  ~ 1 .3 /(C5/C2)  \fC  has  been  used.  This  should  be  considered 
in  the  design  of  the  structures,  as  resonance  response  can  prove  to  be  catastrophic.  An  example 
of  the  use  of  this  method  would  be  the  following: 


f = Sx106  kN/m2 
uc  = 1 .0  m/s 
a = 45° 

u = o.t 


p = 0.93  t/m3 
e/d  = 0.1 
13  = 60° 
e = 0.2 


C = 0.16 


C = 0.16 


J~jO  — ! — £l  a 

f pul  sina  C2 


1 


5 ■ 106  

0.93  1.02  sina  c2 


C?  = 371  -4-  -J-  C? 
J sina  C*  J 


c.  = 1 - p =1-0.1  = 0.755 

1 sina  0.707 

C,  =p+iM  = o.1  + 1=231=2.549 

1 sina  0.707 


C, 

=0.296 

c2 


= c,  = 6 + 6^  cosa  = 1 .777  + 0.424  = 2.201 

C2  ■*  C2  d 

C = 371  • 0.296  • 2.22  = 751 .8. 

Engineering  formula  (eq  33) 

CF=i*l#-  = 0.n. 

x/c  V/75G8 
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r iqure  1 5.  Comparison  with  Korzhavin 's  results. 


COMPARISON  WITH  KORZHAVIN’S  RESULTS 

In  Korzhavin  (1971,  Table  34,  p.  165)  he  has  given  the50  values  for  the  reduction  coefficients. 
In  order  to  compare  the  S0  values  with  results  given  in  this  paper,  they  have  been  multiplied  with 
tan(3,  using  formula  1 1 .30  in  Korzhavin’s  report. 

The  results  have  been  plotted  in  Figure  1 5. 

Detailed  information  concerning  the  system  constants  used  by  Korzhavin  is  not  available.  The 
following  constants  have  been  used  in  the  curve  for  /3  = 60°:  t = 5 x 1 06  kN/m2,  p = 0.93  tonne/m3, 
ue  = 1 .5  m/s,  e/d  0.1 , p = 0.1 , and  e = 0.2. 

The  reduction  coefficient  CF  has  been  transformed  to  CF  K = CF/(2CF  + 0.2)  as  Korzhavin  is 
using  the  bending  strength  in  his  formula.  As  can  be  seen  in  Figure  1 5,  the  plotted  curve  agrees 
fairly  well  with  the  results  ol  Korzhavin. 


RECOMMENDATIONS  FOR  FUTURE  RESEARCH 

In  a work  of  this  kind  there  will  always  be  some  unsolved  problems  which  require  further  study. 

It  has  been  possible  to  develop  a theory  whereby  the  forces  can  be  predicted  with  an  acceptable 
accuracy,  and  model  tests  have  been  made  in  order  to  verify  the  theory.  However,  certain  problems 
are  not  solved,  and  therefore  it  has  been  necessary  to  use  empirical  solutions  in  some  cases.  The 
width  of  the  rupture  channel  has  not  been  theoretically  determined,  and  further  developments  should 
make  this  possible.  The  combined  axial  and  bending  failure  also  needs  further  investigation.  The 
rto/rc  ratio  should  be  reevaluated  as  this  is  quite  dependent  on  the  method  of  performing  compression 
tests  (the  IAHR  ice  committee  is  at  present  studying  this  problem). 

The  hydrodynamic  mass  of  the  ice  debris  should  be  further  evaluated  by  model  tests.  The 
rupture  pattern  and  the  irregularity  of  the  rupture  lines  should  also  be  examined  as  this  has  an 
influence  on  the  dynamic  forces.  Finally,  it  is  emphasized  that  model  tests  with  artificial  as  well  as 
natural  ice  should  be  continued  in  order  to  clarify  the  remaining  unsolved  problems  in  this  paper. 
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APPENDIX  A.  TEST  DATA 
Introduction 

Some  of  the  test  results  are  given  in  the  following  tables.  The  bending  strength  rb  and  the 
compression  strength  rc  of  the  ice  varied  from  test  to  test,  but  in  this  presentation  the  strength 
parameters  have  all  been  normalized  to  a mean  value  using  e = 0.22.  All  the  force  registrations 
have  been  corrected  accordingly.  The  rupture  distance y is  determined  from  y = x sina,  where  x 
is  measured  in  the  direction  of  the  motion. 

Detailed  information  from  the  test  is  available  in  an  ISVA  internal  report  by  C.  Sorensen  (1974). 
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Te«  with  artificial  ice:  Mo.  3/7.1  Test  with  artificial  ice:  No.  3/7. 
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Test  with  artificial  ice:  No.  3/8.1  Test  with  artificial  lc..  No.  3/9 
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Test  with  artificial  lea:  No.  3/9.2  Test  with  artificial  ice:  No.  3/10. 
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